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Abstract.  We  present  analytic  models  for  the  steady  state 
potential  distributions  surrounding  a  spinning,  dielectric- 
coated,  spherical  spacecraft  charging  in  sunlight.  The  sun 
direction  is  assumed  to  lie  in  the  satellite  bellyband  plane, 
perpendicular  to  the  spin  axes.  The  models  are  based  on  a 
multipole  expansion  of  Laplacian  potentials  external  to  the 
spacecraft  surface.  The  combination  of  monopole  potentials 
along  with  the  dipole  or  quadrupole  contributions  produce 
potential  barriers  which  form  at  the  satellite  surface.  These 
barriers  can  block  escaping  photoelectrons  and  lead  to  cur¬ 
rent  balance,  allowing  sunlight  charging  to  high  negative  lev¬ 
els.  In  a  previous  treatment,  analytic  models  were  limited  to 
fast  spin  relative  to  differential  charging  rates  so  that  the  so¬ 
lutions  had  azimuthal  symmetry  around  the  spin  axes.  By  in¬ 
troducing  an  associated  Legendre  term  into  the  potential  ex¬ 
pansion,  the  azimuthal  symmetry  is  removed,  and  the  models 
can  be  developed  to  encompass  any  spin  rate.  The  analysis 
turns  up  three  functions  of  spin  rate  which  are  only  known 
at  the  spin  limits,  but  the  characteristics  of  the  charging  of 
a  rotating  sphere  can  be  explored  using  approximate  forms 
which  represent  the  basic  trends.  For  finite  spin,  the  sunlit 
side  charges  less  (negatively)  than  the  shade  side  which  is 
in  contrast  to  the  fast  spin  case,  where  these  two  potentials 
are  equal.  Also,  for  finite  spin,  differential  charging  devel¬ 
ops  perpendicular  to  the  sun  and  spin  axis  directions,  due 
to  the  transverse  motion.  This  transverse  charging  occurs  at 
all  finite  spin  rates,  disappearing  only  at  the  zero  and  infi¬ 
nite  spin  limits.  There  is  a  correlated  lag  angle  between  the 
direction  of  maximum  sheath  radius  and  the  sun  line.  Plots 
are  given  to  illustrate  the  potential  distributions  representing 
barrier  dominated  sunlight  charging  of  a  spinning  dielectric 
coated  spherical  satellite. 

Keywords.  Space  plasma  physics  (Electrostatic  structures; 
spacecraft  sheath,  wakes,  charging;  General  or  miscella¬ 
neous) 
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1  Introduction 

It  is  known  that  a  satellite  in  sunlight  can  charge  to  substan¬ 
tial  negative  voltages,  even  though  the  photoelectron  current 
(positive)  from  surfaces  (Hinteregger,  et  al.,  1959;  Wrenn 
and  Heikkila,  1973)  is  typically  much  larger  than  the  ambient 
currents  in  space.  This  effect  occurs  because  photoelectrons 
escaping  from  the  satellite  can  be  blocked  by  potential  barri¬ 
ers  at  the  surface,  leading  to  current  balance  and  to  charging 
to  high  (negative)  values.  In  this  paper,  we  consider  barrier 
dominated  charging  models  for  spinning  spherical  satellites. 
Only  steady  state  configurations  are  treated,  although  we  use 
a  time  dynamic  argument  to  obtain  time  scales  for  the  charg¬ 
ing  behavior  during  a  spin  period. 

The  analytic  models  assume  Laplacian  potentials  exterior 
to  the  spherical  satellite  surface.  The  potentials  are  expanded 
in  a  series  of  multipoles  and  we  consider  the  three  lowest  or¬ 
der  terms:  the  monopole,  dipole,  and  the  quadrupole.  The 
monopole-dipole  model,  where  the  satellite  is  non-rotating, 
has  been  treated  by  Besse  and  Rubin  (1980),  Mandell  et 
al.  (1978),  and  Higgins  (1978).  The  monopole-quadrupole 
system  has  been  treated  by  Tautz  (2003).  The  monopole- 
dipole  model  forms  with  the  dipole  aligned  with  the  sun  di¬ 
rection.  In  the  monopole-quadrupole  system,  the  sun  line 
lies  in  the  bellyband  plane  and  the  satellite  is  spinning  fast. 
A  combination  of  these  models,  representing  charging  with 
an  arbitrary  sun  direction  with  respect  to  the  spin  pole,  was 
developed  by  Tautz  and  Lai  (2005).  In  that  treatment  approx¬ 
imate  azimuthal  symmetry  was  assumed,  which  would  be  the 
case  if  the  satellite  is  rapidly  rotating.  Here  we  drop  this  re¬ 
striction  and  allow  the  spin  rate  to  vary  from  zero  to  the  fast 
spin  limit,  which  enables  the  models  to  encompass  a  much 
wider  class  of  satellites.  To  simplify  the  analysis,  we  assume 
that  the  sun  direction  lies  in  the  bellyband  plane. 

The  basic  equations  describing  the  models  and  their  pa¬ 
rameters  are  given  in  Sect.  2.  A  criterion  for  for  slow  or  fast 
spin  is  discussed  in  Sect.  3.  Section  4  gives  the  calculation  of 
the  photoemission  barrier.  In  Sect.  5,  approximate  solutions 
for  a  variable  finite  spin  rate  are  outlined  and  in  Sect.  6  plots 
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are  shown  to  describe  the  characteristic  features  of  a  spinning 
system.  Section  7  contains  a  summary  and  conclusions. 

2  Description  of  the  models 

Consider  a  dielectric  covered  spherical  satellite  that  is  spin¬ 
ning  in  sunlight.  If  the  ambient  charge  density  is  low,  which 
is  a  typical  condition  at  geosynchronous  altitudes,  the  po¬ 
tentials  outside  such  a  satellite  are  given  approximately  by 
a  solution  to  Laplace’s  equation.  In  spherical  coordinates, 
the  exterior  Laplacian  potentials  are  of  the  form  (see  Magid, 
1972): 

tt 

n 

E  cm  n  m  cos(wi^)  +  Cnm  sm(m(p)]  I  (1) 

m  J 

where  r  is  the  radius,  0  is  the  polar  angle  and  (p  is  the  azimuth 
angle.  The  P™  ( x )  are  the  associated  Legendre  polynomials 
(Beyer,  1973)  with  x=cos(0).  The  sum  is  over  n=0,  1, 
2. . .  and  the  sum  £m  is  over  m=0  to  n.  The  constant  coeffi¬ 
cients  Bnm,  Cnm  depend  on  the  exterior  conditions. 

In  terms  of  Cartesian  coordinates  (X,Y,Z)  in  the  space 
frame  of  reference,  the  satellite  is  assumed  to  be  spinning 
around  the  Z-axis  with  the  sun  direction  lying  along  the  pos¬ 
itive  X-axis.  In  this  configuration,  the  potentials  should,  by 
symmetry,  be  even  in  x  (i.e.  with  respect  to  +Z  or  — Z).  We 
can  eliminate  the  Bi\  and  C21  terms  because  P\(x)  's  odd  in 
x  (the  Bio  term  could  also  be  omitted  at  this  point,  but  we 
leave  it  in  to  show  the  form  of  the  equations).  The  B22  term 
can  be  dropped  because,  when  it  is  expressed  in  Cartesian 
coordinates,  it  is  even  in  the  X-coordinate  and  we  expect  to 
see  sunlit  to  shade  side  differences  for  the  sun  on  one  side 
at  +X.  The  C22  term  can  also  be  dropped  for  the  same  rea¬ 
son,  since  it  is  invariant  under  X,Y  reflections.  The  potential 
expansion  for  n=0,  1 ,  2  then  has  the  form 

V(r,B,<f>)  =  j  {b0°(x)+  i  [pfaMio  +  />,'(*) 

(A  1 1  cos (p  +  A sin 0)]  +  -^P2°(x)A2o  +  ...  j  (2) 

Here  the  first  term  is  the  monopole,  the  next  three  are 
Z,  X  and  Y  dipoles,  and  the  last  is  the  quadrupole  com¬ 
ponent.  All  higher  order  terms  are  neglected.  For  con¬ 
venience,  in  Eq.  (2),  we  have  factored  out  the  first  term 
so  that  K=Boo  represents  the  monopole  potential  and  the 
coefficientsA„m=Bm„/X,  Anm*=Cnm/K  give  the  strength 
relative  to  the  monopole.  In  this  expression  the  Legendre 
polynomials  are  (Beyer,  1973  ) 

C  =  1 ;  P?  (-*)  =  -v ;  p,1  (x)  =  (1  -  -v2) ' /2 ; 

P2°(x)  =  l-  (3.x2  -  1)  (3) 


The  Pj1  function  was  absent  from  previous  treatments  be¬ 
cause  the  potential  was  assumed  to  have  azimuthal  symmetry 
(m=0)  and  it  is  this  P,1  term  which  enables  consideration  of 
the  finite  spin  cases. 

The  K  parameter  in  the  models  depends  on  the  balance 
of  the  incoming  and  outgoing  satellite  surface  currents  and 
weakly  on  the  photoemission  barrier  region  (photosheath). 
For  negative  charging,  K  is  less  than  zero.  The  Anm,  Anm* 
parameters,  which  give  the  strength  of  the  non-monopole 
contributions  to  the  expansion,  are  also  used  in  setting  cur¬ 
rent  balance  and  they  depend  more  sensitively  on  the  photo¬ 
sheath  structure.  To  actually  determine  these  parameters  one 
needs  to  know  the  environment  currents,  the  satellite  surface 
material  properties  and  the  photoemission  model.  These  in¬ 
puts,  which  can  be  quite  complicated,  are  not  considered  in 
our  schematic  approach  and  we  take  K  and  Anm ,  Anm  *  to  be 
free  parameters  of  the  models. 

If  we  normalize  to  the  sphere  radius  so  that  at  the  satellite 
surface,  r= 1,  and  substitute  for  the  Legendre  terms,  Eq.  (2) 
becomes: 

V(\ ,  0,  (p)  =  K  {l  +  cosP  Aio+sin0  (An  cos^+A,  *  sin <p) 
+  I(3cos20-i)a2o+  }  (4) 

For  a  satellite  spinning  around  the  Z-axis  with  the  sun  di¬ 
rection  towards  positive  X,  the  surface  potentials  at  the  co¬ 
ordinate  nodes  +Z,  -Z,  +X,  -X,  +Y,  -Y  can  be  labeled  as 


follows: 

VN  =  V(l,  0°,  any)  =  K  (1  +  A,0  +  A20) 

(north  spin  pole)  (5) 

=  V(\,  180°,  any)  =  K  (1  -  A10  +  A20) 

(south  spin  pole)  (6) 

VF  =  V  (1, 90°,  0°)  =  K  (1  +  An  -  A20/2) 

(front,  sunlit  side)  (7) 

VB  =  V  (1,  90°,  180°)  =  K  (1  -  AU  -  A20/2) 

(back,  shaded  side)  (8) 

VE  =  Y(1,900,900)  =  X(1  +Au*-A2o/2) 

(east  side)  (9) 

VW  =  V(l,  90°,  -90°)  =  K(  1  -  An*  -  A20/ 2) 

(west  side)  (10) 

From  the  sum  and  difference  of  Eqs.  (5),  (6)  we  get 

(Vs  -  Vn)/2  =  K(Aio)  (11) 

(Vs  +  VN)/2  =  K(\  +  A20)  =  VP 

(average  pole  potential)  (12) 
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It  is  apparent  that  if  Aio=0,  then  Vs—Vn-  By  symmetry,  this 
will  only  occur  when  the  sun  lies  in  the  X-Y  plane  (belly- 
band).  In  this  paper,  the  sun  direction  is  assumed  to  lie  along 
the  X-axis  i.e.  we  will  fix  Aio=0.  An  approximate  way  to 
treat  more  general  sun  directions  in  the  fast  spin  limit  has 
been  given  in  Tautz  and  Lai  (2005). 

We  now  take  the  sum  and  difference  of  Eqs.  (7),  (8)  and 


Eqs.  (9),  (10)  to  get 

(VF-VB)/2  =  KAu=-D/2  (13) 

(VF  +  VB)/2  =  K  (1  -  A20/ 2)  (14) 

(VE-  Vw)/2  =  KAl{  =  -d/2  (15) 

(Ve  +  Vw)/2  =  K(\-A2o/2)  (16) 


The  variables,  D  and  d,  denote  the  sun  aligned  and  transverse 
differential  charging  respectively,  and  are  seen  to  be  simply 
related  to  the  X  and  Y  dipole  strengths.  In  these  equations, 
the  signs  of  the  potentials,  including  K,  would  be  negative 
for  typical  charging.  The  potentials,  normalized  to  K,  would 
therefore  be  positive.  For  the  sun  direction  at  +X,  we  ex¬ 
pect  VB  /  K>Ve/K  and  An<0.  For  satellite  spin  in  the 
counter  clockwise  sense,  we  would  have  V\v/K>VE/K,  giv¬ 
ing  Au*<0.  Spin  in  the  clockwise  sense  would  correspond 
toAu*>0. 

By  comparing  Eqs.  (14),  (16)  we  get 

VE  +  Vw  =  VB  +  Vp  (17) 

This  condition  states  that,  for  the  configuration  and  approx¬ 
imation  level  that  we  have  considered,  the  system  is  diago¬ 
nally  balanced.  We  mean  by  this  that  the  voltage  change  in 
any  spin  quadrant  is  the  same  as  the  change  in  the  diagonally 
opposite  quadrant.  This  can  be  seen  in  Fig.  1,  which  gives  a 
top  view  of  the  satellite  spin  quadrants,  with  the  sun  located 
to  the  right  and  the  rotation  counter  clockwise.  In  the  figure, 
the  symbols  E,  W,  B,  F  label  the  east,  west,  back  and  front 
sides  of  the  sphere.  From  Eq.  (17),  it  can  be  seen  that  the 
voltage  increment  in  the  first  quadrant  (Vg-Vp)  equals  the 
decrement  in  the  third  quadrant  (Vw-Vg).  The  same  type  of 
voltage  balance  holds  for  quadrants  two  and  four. 

Equations  (14),  (16)  lead  to 

A 20  =  2(1-  Vm/K)  =  2  (1  -  (VE  +  W)/2 K) 

=  2  (1  -  (V*  +  VV)/2 K)  (18) 

where  we  denote  by  Vm  the  average  middle  (bellyband)  po¬ 
tential 

VM  =  -(VE  +  VW  +  VF  +  VB)  (19) 

4 

We  can  see  that  the  quadrupole  strength  coefficient,  A?o,  can 
be  determined  by  the  ratio  of  Vm  to  K  and  this  coefficient  is 
zero  if  Vm  =  K. 


E 


Fig.  1.  A  top  view  of  the  satellite  spin  quadrants,  with  the  sun  lo¬ 
cated  to  the  rieht  and  the  rotation  counter  clockwise. 


The  average  surface  potential  is  denoted  by 

<  V  >=-{VN  +  Vs  +  4Vm)  (20) 

0 

From  the  sum  of  Eqs.  (5-10)  we  find  that  <V>—K.  This 
simple  equation  is  helpful  in  understanding  the  overall  charg¬ 
ing  behavior  of  spinning  spheres.  It  can  be  written  in  the 
equivelant  form 

VP/K  =  3-2VM/K  (21) 

We  see  that,  if  the  bellyband  potentials  change,  the  spin  pole 
potentials  adjust  linearly.  In  general,  any  change  to  a  surface 
potential  will  be  compensated  by  changes  to  other  surfaces 
such  that  the  average  surface  potential  remains  at  K. 

If  we  are  given  the  six  node  potentials,  the  coefficents  Aio, 
An,  An*  and  A2o  are  determined  by  Eqs.  (11),  (13),  (15), 
(18).  We  now  examine  two  limit  cases  for  which  these  node 
potentials  are  known.  We  first  look  at  the  zero  spin  limit.  It  is 
known  that  the  solution  is  a  monopole-dipole  system  (Besse 
and  Rubin,  1980)  pointing  in  the  sun  direction,  X ,  and  given 
by 

VB  =  K(l  +A) 

VF  =  K{  1  -  A) 

VN  =  Vs  =  VE=  Vw  =  K  (22) 

where  A>  '/2  is  the  dipole  strength  parameter.  Solving  for  the 
coefficients,  we  find  Aio=An*=A2o=0  and  An=— A.  Here 
Vm=K  and  this  makes  the  quadrupole  term  go  to  zero.  The 
differential  charging  is  given  by  D=2KA  and  d= 0.  There  is 
also  differential  charging  between  the  spin  poles  and  belly¬ 
band  which  is  azimuth  dependent.  We  next  look  at  the  case 
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Fig.  2.  Hypothetical  linear  charging  voltages  as  functions  of  the 
spin  period  fractional  times.  The  charging  rate  on  the  dark  side  is 
assumed  to  be  equal  to  the  discharge  rate  on  the  sunlit  side.  The 
horizontal  dashed  line  represents  the  dark  side  equilibrium  value. 
The  top  letters  denote  the  east  (E),  back  (B),  west  (W),  and  front 
(F)  surfaces  relative  to  the  sun  direction. 

of  infinite  spin.  It  is  known  that  the  solution  is  a  monopole- 
quadrupole  system  (Tautz  and  Lai,  2005)  given  by 

VE=Vw=VF=VB=K{\-A/2) 

Vn=Vs=K(.\+A)  (23) 

where  A  >2/3  is  now  the  quadrupole  strength  parameter. 
Again,  solving  for  the  coefficents  we  get  A  io= A 1 1 = A  n  *=0 
and  A20—A.  The  differential  charging  is  here  D=d~ 0  and 
all  the  bellyband  potentials  are  equal,  which  characterizes 
the  fast  spin  limit.  There  is  again  pole-bellyband  differential 
charging,  which  in  this  case  is  azimuthally  symmetric. 

We  have  seen  that,  depending  on  the  parameters  A  io,  A  n , 
A  n  *  and  A20,  we  can  get  a  zero  spin  or  infinite  spin  solution. 
In  both  cases  A  io=0  which  signifies  that  the  sun  lies  in  the  x- 
y-plane.  Also,  in  both  cases,  we  have  Au*=0,  but  we  do  not 
assume  that  this  is  true  in  general. 

3  Fast  and  slow  spin 

We  now  wish  to  discuss  the  models  at  finite  spin  rates,  be¬ 
tween  the  limit  cases  of  zero  and  infinity.  To  get  an  insight 
into  the  charging  time  scales,  we  need  a  means  of  judging 
whether  a  given  spin  rate  is  fast  or  slow.  We  can  use  a  time 
dynamics  argument  to  characterize  the  significant  charging 
intervals  during  a  rotation.  The  analysis  takes  into  account 
that  each  surface  element  can  be  characterized  by  its  capaci¬ 
tance  per  unit  area  and  it  will  charge  and  discharge  on  a  time 
scale  set  by  that  value  as  well  as  by  the  spin  rate.  There  are 
thus  two  relevant  times  to  compare:  the  time  for  equilibrium 
surface  charging,  T,  and  the  satellite  spin  period,  r.  We  can 
define  a  spin  rate  parameter,  R,  as 

R  =  2T /r  (24) 

The  factor  of  two  is  inserted  because  it  is  the  spin  half  pe¬ 
riod  which  is  the  critical  time.  It  should  be  noted  that  R  is 


only  defined  approximately.  It  depends  on  the  charging  time 
T,  which  is  not  precisely  defined.  T  can  be  estimated  from 
the  time  constant  for  differential  charging  ,  but  there  is  not  a 
precise  value  at  which  a  transition  to  the  fast  spin  state  is  real¬ 
ized.  Experimentally,  the  fast  spin  limit  could  be  recognized 
by  a  lack  of  spin  modulation  in  measured  on-board  surface 
fluxes. 

Consider  the  behavior  of  an  isolated  surface  cell  on  a 
spinning  spacecraft.  If  the  satellite  is  in  a  charging  en¬ 
vironment  (Lai  and  Tautz,  2006)  a  surface  on  the  shaded 
side  will  go  negative  at  a  rate  dWd t=J/C  where  J  is  the 
dark  side  current  density  and  C  is  the  capacitance  per  unit 
area  of  the  element.  Typical  values  for  an  orbiting  geosyn¬ 
chronous  spacecraft  are  7=1  jiAJm2  and  C=1  /zF/m2  which 
gives  dV7d/=l  Volt/s.  If  the  satellite  spins  at  one  RPM,  the 
charging  on  the  dark  side  would  reach  30  Volts.  This  is  not  a 
negligible  amount  for  a  spacecraft  charging  to  a  few  hundred 
volts  in  sunlight. 

To  estimate  time  scales,  we  follow  a  cell  as  it  moves 
around  the  spacecraft  bellyband  during  one  spin  period 
through  surface  points  E,  B,  W,  F  and  then  back  to  E  (see 
Fig.  1).  The  corresponding  fractional  times  are  0,  1/4,  1/2, 
3/4  and  1.  In  Fig.  2,  we  show  five  hypothetical  charging 
curves  as  a  function  of  the  spin  period  fractional  times.  For 
simplicity,  the  curves  are  shown  as  linear.  The  charging  rate 
on  the  dark  side  is  assumed  to  be  equal  to  the  discharge 
rate  on  the  sunlit  side.  The  curves  are  drawn  as  diagonally 
balanced,  in  the  sense  discussed  in  Sect.  2.  The  horizontal 
dashed  line  represents  the  dark  side  equilibrium  level,  Ve.  In 
the  figure,  the  curves  are  labeled  1-5  and  locate  the  main 
regions: 

1 )  R= 0  -  the  stationary  case 

2)  R=\I2  -  the  end  of  the  very  slow  spin  region 

3)  /?=1  -  the  end  of  the  slow  spin  region  and  the  start  of 
the  fast  spin  region 

4)  R= 2  -  the  center  of  the  fast  spin  region 

5)  /?= 00  -  a  horizontal  line,  which  is  the  fast  spin  limit 

The  parameter  value  /?=1  specifies  the  boundary  between 
slow  and  fast  spin.  The  charging  behavior  is  quite  differ¬ 
ent  between  these  two  regions.  For  slow  spin,  the  potential 
levels  off  at  the  dark  side  equilibrium  level,  Ve.  For  fast  spin, 
the  dark  side  potential  never  gets  to  equilbium.  It  is  cut  off  at 
time  r/2,  at  a  value  which  decreases  with  increasing  spin. 

As  in  Sect.  2,  we  denote  the  front  to  back  voltage  differ¬ 
ence  in  the  sun  direction  as  D  and  the  difference  in  the  trans¬ 
verse  direction  as  d.  From  Fig.  2,  we  can  see  that  the  differ¬ 
ential  charging,  D,  is  a  maximum  at  zero  spin  and  decreases 
with  R.  At  7?=1,  the  average  potential  on  the  dark  side  ap¬ 
proximately  equals  the  average  on  the  discharge  side,  and  D 
goes  to  zero,  where  it  remains  for  R>\.  The  transverse  dif¬ 
ferential  charging,  d.  on  the  other  hand,  occurs  whenever  the 
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satellite  has  finite  spin.  It  is  due  to  the  charging  up  of  the 
cell  as  it  moves  from  east  to  west  and  then  discharges  from 
west  to  east,  creating  a  transverse  voltage  difference.  It  only 
vanishes  at  the  limits  when  the  satellite  is  stationary,  so  that 
by  symmetry  the  potentials  are  the  same,  or  at  the  fast  spin 
limit  where  all  bellyband  potentials  are  equal. 

Real  charging  curves  would  be  non-linear.  Any  curvature 
would  smooth  out  the  linear  results.  For  example,  hard  zeros 
for  D  would  go  into  small  numbers.  The  time  dynamics  anal¬ 
ysis  is  schematic,  as  the  models  pertain  solely  to  the  steady 
state  and  the  surface  is  not  composed  of  isolated  cells,  but  it 
serves  to  identify  the  charging  regions. 

To  illustrate  the  effect  of  time  averaging,  we  calculate  the 
potentials  in  the  circulating  cell  case,  which  can  be  easily  ob¬ 
tained  for  the  linear  curves  of  Fig.  2.  The  results  are  shown 
in  Table  1 .  Here  the  time  averages  are  taken  using  plus  or 
minus  one  quarter  of  a  spin  period,  centered  on  the  fractional 
time  for  the  potential.  For  this  calculation,  we  have  used  po¬ 
tential  limits  v^=3/2  and  U£=l/2,  which  for  comparison  have 
been  chosen  to  give  the  Besse-Rubin  threshold  potential  ra¬ 
tios,  V/K,  at  zero  spin  (see  line  1  of  Table  1).  The  computed 
time  averaged  charging  variables  are  summarised  in  the  ta¬ 
ble.  The  charging  behavior  is  similar  to  the  dynamic  cases. 
We  note  that  the  time  averaged  potentials  Ve,  Vw  are  equal 
at  R= 0,  rather  than  having  discontinuous  jumps. 

We  can  parameterize  the  spin  behavior,  by  letting  the 
model  coefficients  An,  A1p  A20,  be  functions  of  R.  We 
will  assume  these  functions  have  the  form 

Anm(R)  =aanm(R)  (25) 

and  a  similar  form  for  A 1 1  *.  Here  or  is  an  overall  strength  pa¬ 
rameter  and  the  anm(R  )  functions  give  the  relative  weights. 
If  a  is  zero,  there  is  no  differential  charging  and  we  revert  to 
the  eclipse  charging  case.  Using  the  K,  a,  R  parameteriza¬ 
tion,  the  expression  for  V  becomes 

V(r,0,4>)  =  — 

r 

(ot 

1  H — sin#  [au(/?)cos$  +  ax*(R)  sirup] 

+  ^2  (3cos2#  —  l)a20(fl)}  (26) 

This  is  the  approximate  potential  that  is  used  to  describe  a 
spinning  sphere.  The  A 10  term  has  been  set  to  zero.  The  R 
dependence  is  discussed  more  fully  in  the  following  sections. 


4  The  photoemission  barrier 


Table  1.  Time  averaged  potentials  for  a  rotating  surface  cell. 
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vB 

VF 

% 

VE 

D 

d 

0 

3/2 

1/2 

1 

1 

1 

0 

1/4 

11/8 

5/8 

9/8 

7/8 

3/4 

1/4 

1/2 

5/4 

3/4 

5/4 

3/4 

1/2 

1/2 

3/4 

9/8 

7/8 

31/24 

17/24 

1/4 

7/12 

1 

1 

1 

5/4 

3/4 

0 

1/2 

3/2 

5/6 

5/6 

1 

2/3 

0 

1/3 

2 

3/4 

3/4 

7/8 

5/8 

0 

1/4 

4 

5/8 

5/8 

11/16 

9/16 

0 

1/8 

OO 

1/2 

1/2 

1/2 

1/2 

0 

0 

has  a  solution.  We  can  write  the  potential,  Eq.  (26),  in  sim¬ 
plified  form  as 


V(r,  9,  <p)  = 


K 


aT\  uTi 
1  +  —  +  -==■ 
r  rl 


where  we  have  defined  the  terms 
Tj^sin#  [aii(^)cos0  +  a,*(R)sin</>] 

(3  cos 2  9  -  l) 


72  = 


L«20  (R) 


(28) 


(29) 

(30) 


We  get  from  Eqs.  (27),  (28)  a  quadratic  equation  for  r 
r2  +  2aTi  r  +  3aT2  =  0  (31) 


The  solution  rb  for  the  barrier  radius  is  thus 


r  ■>  *i  1/2 

rb  =  -«r,+  [(crW,)2  -3«r2J  (32) 

Note  that  rb  depends  on  a  and  0,  <p,  R  (through  T\  and  73), 
but  not  on  K.  Since  the  K  parameter  does  not  appear  in 
Eq.  (32),  it  is  mainly  the  a  parameter  which  represents  the 
photosheath  structure. 

We  can  write  rb  =  1  +  s,  where  s  is  the  sheath  width,  nor¬ 
malized  to  the  sphere  radius.  If  the  sheath  forms  close  to 
the  sphere  surface  (^<^1),  we  find  that,  neglecting  quadratic 
terms,  we  have 


1  +2g:7’i  +3a72 
2  (1  +  aT\) 


(33) 


The  condition  giving  the  threshold  for  charging  (rb=  1,  5=0) 
is 


1  +  20:7)  +  30:73  =  0 


(34) 


For  substantial  negative  charging  to  occur  in  sunlight  it  is  ex¬ 
pected  that  a  potential  barrier  will  form  just  outside  the  sur¬ 
face  and  trap  escaping  photoelectrons.  Such  a  barrier  exists 
if 


d_ 

Tr 


v(r,e,4>)  =  0 


(27) 


Taking  a=l  at  the  threshold,  we  find  that  the  values  of  T 1 ,  T2 
in  the  zero  spin  case  are  T|=— 1/2,  T2=0  and  in  the  fast  spin 
limit  are  Tj=0,  T2=— 1/3. 

For  a  non-spinning  spacecraft  we  expect  that  the  max¬ 
imum  sheath  radius  would  occur  along  the  sun  direction 
(#=90°,  <p= 0°).  However,  for  a  spinning  satellite  this  is  not 
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necessarily  true.  To  get  the  location  of  the  maximum  sheath 
radius,  we  differentiate  Eq.  (32)  with  respect  to  <p.  Since  the 
4>  dependence  occurs  only  through  T\,  we  can  use  dT\/d(j>=0, 
which  leads  to  the  sheath  angle,  <pb,  at  the  maximum  barrier 
radius,  as 

^tan-'^W/nuC/?)]  05) 

If  we  evaluate  the  threshold  condition  along  the  line  of  max¬ 
imum  sheath  radius  (0=90°,  <p=<pb),  we  get  from  Eqs.  (29), 
(30),  (34),  (35)  a  relation  between  the  coefficients 

3  1 

-  an(R)  cos  (pb  —  a{*(R)  sin<pb  +  -a2o(R)  =  -  (36) 

4  2 

As  R  goes  from  0  to  co,  the  coefficients  will  change,  but 
they  will  always  satisify  Eq.  (36)  at  threshold.  For  or>l, 
the  Anm(R)  coefficients  would  be  scaled  from  the  anm(R ) 
values,  resulting  in  a  finite  barrier  radius  (rb>  1)  outside  the 
surface,  with  possible  charging. 

The  model  barrier  height,  Vb,  is  defined  as 


Vb  =  V(r„,9^)-V{ \,e,<t>) 


which  yields 


jl-r*  +  a7i 


0-# 

rb 


+  aT2 


(37) 


In  general,  the  ratio  Vb/K  depends  on  a  and  0,  <p,  R.  At 
the  maximum  sheath  radius,  7)  should  be  evaluated  using 

4>=<Pb(R  )• 

At  R= 0,  Eqs.  (32),  (37)  give  the  characteristic  equations 
for  a  monopole-dipole  system  (Besse  and  Rubin,  1980) 


rb  =  a 


Vb  =  ^(a-\)2  (38) 

2a 

and  as  R-+co  they  go  into  the  fast  spin  equations  (Tautz  and 
Lai,  2005) 


Vb  =  K  (ja  +  joT1'2  -  l)  (39) 

Note  that  these  expressions  differ  from  ones  given  in  former 
publications  because  the  definition  of  a  has  changed  slightly 
from  the  previously  used  non-monopole  strength  parameter. 
We  can  see  directly  that  a=l  gives  back  the  threshold  values 
rb=  1,  Vb/K=0. 

The  energy  of  the  emitted  photoelectrons  is  low  (1  to  2  eV, 
Whipple,  1981;  Lai  et  al.,  1986),  and  the  barrier  height  Vb 
required  to  stop  a  substantial  fraction  of  them  would  only 
be  a  few  volts.  For  high-level  charging,  with  K  at  hundreds 
of  volts,  the  normalized  barrier  height  Vb/K  will  be  near  to 
its  threshold  value.  It  is  this  approximate  feature  of  sunlight 


charging  which  allows  the  threshold  states  to  be  used  effec¬ 
tively  as  a  basis  for  the  solutions. 

The  above  analysis  of  the  photosheath  is  based  wholely 
on  Laplace’s  equation.  To  do  the  photoelectron  dynamics 
consistently  would  require  solving  Poisson’s  equation,  with 
calculated  electron  density  in  the  sheath  region,  which  is  be¬ 
yond  the  scope  of  our  schematic  models.  A  photosheath 
model  could  be  added  separately  and  would  provide  a  re¬ 
lation  between  the  barrier  height,  Vb,  and  the  sheath  width 
s.  This,  combined  with  the  current  balance  condition,  would 
lead  to  a  consistant  closed  system  (see  the  comment  at  the 
end  of  Sect.  5).  However,  it  can  be  shown  that  for  strong 
differential  charging,  the  space  charge  in  the  photosheath  has 
only  a  small  effect  on  the  barrier  (Mandell  et  al.,  1978)  and 
Laplace’s  equation  is  sufficient. 


5  Approximate  solutions  to  the  models 

In  order  to  explore  the  charging  behavior  of  a  satellite  rotat¬ 
ing  at  finite  spin,  we  assume  approximate  functional  forms 
for  the  spin  dependent  threshold  coefficents  a  1 1(7?),  zzi  i  *  (7?), 
a2o(R).  The  functions  are  estimated  as  follows: 

1 )  We  can  get  an  idea  as  to  the  approximate  behavior  of 
the  «ii(R  )  function  by  examining  Fig.  2  and  Table  1.  It  is 
apparent  that  at  R= 0,  (curve  1  in  the  figure)  that  the  differ¬ 
ential  charging,  D,  has  a  maximum  value.  At  R=l,  (curve  3) 
the  value  of  D  has  dropped  down  to  zero  because  the  front 
and  back  potentials  are  equal  and  this  condition  remains  for 
all  fast  spin  rates  (curves  >3).  We  expect  the  dipole  response 
of  a  rotating  sphere  to  be  similar,  but  somewhat  smoother. 
An  approximate  function  that  ramps  down  from  the  X-dipole 
solution  is  D=Xexp(-b/f)  and  from  Eq.  (13) 

a,i(R)  =  - D/(2K )  =  -(1/2)  exp  {-bR)  (40) 

The  factor  of  b  in  the  exponent  is  arbitrary.  A  larger  b  puts 
the  value  of  D  at  R=1  closer  to  zero.  In  the  examples  below, 
we  use  £>=2.0. 

2)  The  behavior  of  the  a\\*(R)  function  can  also  be  esti¬ 
mated  from  Fig.  2  and  Table  1 .  At  R= 0,  the  east  and  west  po¬ 
tentials  are  equal  by  symmetry  (curve  1)  and  d= 0.  At  R=co, 
(curve  5)  we  again  have  d= 0  because  all  bellyband  poten¬ 
tials  are  equal.  For  intermediate  spins,  d  would  go  through 
a  maximum.  We  first  consider  the  region  £?>1.  We  can  see 
that  d  decreases  with  R,  since  it  is  cut  off  at  r/2,  so  that  the 
functional  form  is  y2(R)=c  (x l(2T))=clR  where  c  is  a  con¬ 
stant.  For  the  region  R<\,  the  differential  charging,  d,  starts 
from  zero,  ramps  up,  goes  through  a  maximum  somewhere 
before  R=  1,  and  then  connects  with  the  function  at  /?>  1 . 
We  can  implement  this  behavior  smoothly  with  a  second 
degree  polynomial,  yi(R),  satisfying:  yi(0)=0,  yi(l)=y2(l), 
dy\(\ )/dR=dy2(\ )/dR  which  yields  yi=c  R  (3-2  R).  The  co¬ 
efficient  at  R  <  1  is  then  determined  using  y\  in  Eq.  ( 1 5) 

«,,*(/?)  =  -d/(2K)  =  —c'R{2  —  2R)2)  (41) 
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and  for  R>  1,  using  yi  it  is 


a^(R)  =  -c'/R, 


(42) 


where  c'=c/(2K).  A  value  for  c'  can  be  estimated  by  as¬ 
suming  typical  physical  conditions  at  /J=l.  As  described  in 
Sect.  3,  if  we  have  a  satellite  spinning  at  1  RPM  there  could 
be  about  30  Volts  of  transverse  charging  and,  given  a  rep¬ 
resentative  monopole  charging  level  of  500  Volts,  we  find 
c'=30/1000=0.03. 

3)  Rather  than  introducing  an  independent  estimate  for  the 
aio(R)  function,  we  set  the  quadrupole  coefficient  by  a  phys¬ 
ical  condition.  For  a  sheath  barrier  to  occur  exterior  to  the 
sphere  surface,  the  system  must  be  above  threshold,  speci¬ 
fied  by  condition  Eq.  (36),  which  gives 


a2o(R)  = 


1 


+  au(R)  cos  fo  +a,,*(fl)  sin  <pb 


(43) 


This  expression  has  reasonable  limits  since,  at  R= 0  the 
quadrupole  term  is  equal  to  zero  (see  Eq.  22)  and,  at  R= oo 
it  goes  into  2/3,  the  fast  spin  threshold  value  (as  in  Eq.  23). 

The  above  functions  enable  us  to  determine  potential  co¬ 
efficients  for  a  spinning  sphere  at  the  charging  threshold.  To 
get  the  actual  solution  potentials  one  has  to  scale  these  coeffi¬ 
cients  by  a.  The  functions  represent  smoothly  the  trends  and 
limits.  The  real  functions  are  unknown,  but  these  approxi¬ 
mate  forms  are  sufficient  to  illustrate  the  main  features. 

In  setting  up  these  functions  two  new  “free”  parameters  b 
and  c'  have  been  introduced.  These  parameters  are  not  re¬ 
ally  free,  since  in  a  real  physical  charging  case  they  would 
be  determined  from  material  and  plasma  properties.  They 
are  similar  to  K  and  a  in  that  they  would  be  known  if  the 
problem  was  completely  specified.  For  example,  K  and  a 
could  be  calculated  consistently  as  follows.  The  existance  of 
a  photosheath  model  would  provide  a  relation  B(s)  where  B 
is  the  photoemission  barrier  height.  The  multipole  models 
provide  the  scaling  function  Vb/K=f( s).  Taking  the  ratio  of 
these  two  quantities  would  yield  K(s)=K (a).  The  a  variable 
could  be  varied  to  determine  K  and  thereby  obtain  the  sur¬ 
face  potentials.  Then,  using  a  current  collection  model,  the 
plasma  incident  fluxes  at  the  surface  could  be  calculated.  In¬ 
teractions  of  the  incident  particles  with  the  surface,  using  a 
materials  model,  would  yield  the  net  surface  currents.  The 
value  of  a  could  be  moved  until  current  balance  at  the  satel¬ 
lite  surface  was  obtained  and  this  would  give  a  consistant 
solution. 


6  Characteristics  of  the  models 

To  illustrate  the  R  dependence  of  the  models,  we  show  repre¬ 
sentative  plots.  The  various  limiting  values  for  the  plot  vari¬ 
ables,  at  R= 0  and  R=o o,  are  summarized  in  Table  2.  All 
potentials  shown  are  normalized  to  K  and  are  positive.  In  a 
real  physical  case,  the  unnormalized  potentials  and  K  would 
be  typically  negative. 


Table  2.  The  model  threshold  parameters  at  the  spin  limits. 
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R=o 0 
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Fig.  3.  Representative  plots  showing  the  R  dependence  of  the  abso¬ 
lute  value  of  the  multipole  coefficients.  The  solid  lines  are  for  a=  1 .2 
and  the  dashed  lines  for  alpha=1.0.  The  Aj*  terms  are  multiplied 
by  10  for  clarity.  The  differential  charging  voltages,  D  and  d,  are 
obtained  by  multiplying  An  and  Af  *  by  —  IK. 


In  Fig.  3  we  show  the  assumed  R  dependence  of  the  coef¬ 
ficients  -An,  —An*  (x  10  for  clarity),  A 20.  The  curves  are 
given  for  two  values  of  a,  corresponding  to  threshold  (a=T) 
and  to  a  finite  sheath  radius  (a=1.2).  The  functions  An  and 
A  [  *  are  the  drivers  of  the  models  (A20  is  set  from  them).  The 
amount  of  differential  charging,  D  and  d,  can  be  scaled  from 
the  coefficients  An  and  Au*,  by  multiplying  by  —2 K.  The 
threshold  curves  are  indicated  by  a  dashed  line.  The  term 
—an  ramps  down  from  1/2  to  0.  The  term  — a ,*  ramps  up 
from  0,  goes  through  a  maximum  before  R— 1,  then  relaxes 
back  to  0  at  R= 00.  The  term  020  ramps  up  from  0  and  goes 
to  2/3  at  R=o o. 

In  Fig.  4,  we  show  the  average  bellyband  and  spin  pole 
potentials  at  ar=1.0,  1.2.  The  threshold  curves  are  indicated 
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Fig.  4.  The  average  normalized  bellyband  potential,  V^/K,  and 
the  average  spin  pole  potential,  Vp/K,  versus  R.  The  dashed  lines 
indicate  the  threshold  values. 


0  1  2 

SPIN  RATE  R 


Fig.  5.  The  photoemission  sheath  angle,  </>b,  at  the  maximum  sheath 
radius,  versus  R. 


by  dashed  lines.  This  plot  shows  how  the  potentials  ad¬ 
just  to  changing  spin  rates.  As  the  satellite  spins  faster, 
the  quadrupole  term  increases  from  zero  so  that  V^IK  goes 
lower  via  (from  Eq.  1 8) 


Vm  _  .  _  A.2o(R) 
K  2 


(44) 


According  to  Eq.  (21),  the  average  spin  pole  potential,  Vp, 
must  rise  so  that  the  average  surface  voltage  remains  at  K. 
This  behavior  of  V^/K  and  Vp/K  is  plotted  in  Fig.  4  and 
shows  that  the  average  pole-bellyband  differential  charging 
increases  with  the  spin  rate.  Since  the  sun  aligned  and  trans¬ 
verse  differential  charging  fall  off  with  increasing  spin,  the 
pole-bellyband  charging  is  the  dominant  differential  voltage 
in  the  fast  spin  limit.  The  normalized  difference  between  Vp 
and  Vm  ,  using  Eq.  (21),  is 


(Vp  -  Vm)  (  _Vm\ 

K  \  K) 


(45) 


At  R=  0,  the  difference  is  0,  corresponding  the  X-dipole  solu¬ 
tion,  while  at  R= oo  it  goes  into  the  fast  spin  threshold  value 
of  1 .  This  indicates  that  the  potentials  have  changed  from  a 
monopole-dipole  configuration  to  the  “dumbed”  shape  of  the 
monopole-quadrupole  model.  For  intennediate  spin  rates, 
the  potential  is  formed  from  the  monopole  term,  a  combi¬ 
nation  of  an  X-dipole  and  Y-dipole  which  gives  differential 
charging  in  the  bellyband  plane,  and  the  quadrupole  com¬ 
ponent  which  is  axially  symmetric.  As  R  increases,  the 
quadrupole  term  becomes  stronger  and  brings  in  azimuthal 
smoothing  to  the  solutions. 


In  Fig.  5,  we  show  the  sheath  angle,  cpb,  at  a=1.0  (there 
is  no  variation  with  a).  The  sheath  angle  increases  mono- 
tonically  from  zero  at  R= 0.  At  the  fast  spin  limit  it  goes  to 
90°.  The  slight  wiggle  in  the  curve  is  caused  by  an*  going 
through  a  maximum  near  /?=1. 

Figures  6  and  7  show  the  sheath  radius  and  sheath  height 
at  a=l.l,  1.2  and  1.3.  In  both  of  these  figures,  the  threshold 
corresponds  to  the  abscissa  of  the  plots.  The  figures  indicate 
that,  at  constant  a,  the  sheath  becomes  tighter  and  the  barrier 
height  lessens  as  R  increases.  This  effect  is  most  apparent  at 
R<  1. 

Figures  8  to  1 0  show  the  model  solutions,  for  three  spin 
rates,  in  the  physical  space  surrounding  the  sphere.  The  con¬ 
tour  plots  show  a  T=0  slice  of  data  expressed  in  X,Y,Z  coor¬ 
dinates,  which  are  normalized  to  the  sphere  radius.  In  these 
figures,  the  orientation  symbols  B,  F  would  be  to  the  left  and 
right  and  the  E,  W  symbols  would  be  out  of  the  page.  The  a 
parameter  is  1.2.  The  sun  direction,  as  indicated  by  a  white 
line,  is  to  the  right.  The  plots  are  given  for  R= 0,  l/2,  1 .  In 
the  figures,  one  can  see  the  progression  from  an  X-dipole 
dominated,  zero  spin  case  to  the  quadrupole  dominated  so¬ 
lution  at  fast  spin.  The  most  interesting  case  is  depicted  in 
Fig.  9,  which  shows  the  solution  at  intermediate  spin.  It  is 
evident  that  there  is  a  diminishment  of  the  shaded  side  po¬ 
tentials  relative  to  the  R= 0  solution,  due  to  the  effect  of  az¬ 
imuthal  smoothing  from  the  quadrupole  term. 

Figure  11  depicts  the  R- 1/2  solution  from  above  (Z=0 
plane).  In  this  plot,  the  orientation  symbols  E,  W,  B,  F 
would  be  at  the  same  locations  as  in  Fig.  1.  Since  we  are 
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Fig.  8.  Normalized  potentials  for  the  case  R  =  0.0,  a=1.2,  <pb= 0, 
Y=0. 


Fig.  6.  The  photoemission  sheath  radius/^,  at  (fa,,  versus  R. 


Fig.  7.  The  normalized  photoemission  sheath  barrier  height,  V^/K , 
at  4>b,  versus  R. 


at  intermediate  spin,  there  is  strong  azimuthal  dependence, 
with  relatively  large  front  to  back  differential  charging.  The 
transverse  differential  charging  is  made  evident  by  the  devi¬ 
ation  from  +Y,  —  Y  symmetry.  For  this  case,  the  sheath  angle 
is  4>b=93  degrees,  which  causes  the  shift  in  the  potentials, 
clearly  seen  on  the  forward  side  of  the  sphere.  The  satellite 


VIK 


-3.0  0.0  ^  3.0 


Fig.  9.  Normalized  potentials  for  the  case  R= 0.5,  a=1.2,  $/,= 9.3, 
Y=0. 


spin  is  here  assumed  to  be  counter-clockwise,  but  an  equiv¬ 
alent  view  would  be  that  of  a  stationary  satellite,  with  the 
sun  rotating  around  it  in  the  clockwise  sense.  Looked  at  this 
way,  a  positive  <pb  represents  a  lag  angle  with  respect  to  the 
moving  sun  line. 

We  remark  on  some  limitations  of  the  models  with  re¬ 
gard  to  currents.  In  a  real  charging  case,  satellite  surface 
currents  can  be  quite  complex.  In  the  models,  current  bal¬ 
ance  is  assumed  but  the  spatial  distribution  of  the  surface 
currents  is  not  detailed.  During  charging,  it  is  expected  that 
a  net  negative  current  enters  at  the  shaded  surfaces  (mainly 
electrons,  with  photoemission  absent)  and  a  net  positive 
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Fig.  10.  Normalized  potentials  for  the  case  /?=1.0,  a=  1.2,  </>* =23.9, 
Y=0. 
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Fig.  11.  Normalized  potentials  for  the  case  R= 0.5,  a=1.2,  <pb= 9.3, 
Z=0. 

current  (mainly  the  unsuppressed  high  energy  tail  of  the 
photoemission  spectrum)  occurs  at  the  sunlit  areas.  The  net 
current  at  a  particular  surface  element  would  depend  on  the 
material  interactions  with  the  incident  fluxes  and  these  are 
not  treated  by  the  models.  Another  model  limitation  is  that 
charge  transport  through  the  satellite  by  means  of  electri¬ 
cal  connections  is  not  considered.  Such  conduction  currents 
would  tend  to  reduce  the  degree  of  differential  charging. 

In  the  present  models,  barrier  formation  and  non¬ 
monotonicity  of  potentials  are  due  to  constraints  imposed  by 
the  Legendre  potential  solution  considered.  The  plots  3  to 
1 1  above  represent  charging  configurations  which  are  based 
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on  typical  spacecraft  conditions.  The  models  will  give  a 
better  representation  of  real  physical  charging  as  the  spac- 
eraft/plasma  conditions  approach  those  of  the  theory.  In  or¬ 
der  to  verify  the  theory,  it  is  necessary  to  have  measurements 
of  potential  distributions  on  a  spinning  satellite  with  similar 
conditions  of  the  models.  We  are  not  aware  of  such  measure¬ 
ments.  A  suite  of  computer  simulations  of  sunlight  charging, 
done  over  a  wide  range  of  spins  and  spacecraft/plasma  con¬ 
ditions,  could  be  carried  out  to  verify  the  theory  and  to  better 
establish  the  scope  and  limitations  of  the  models. 

7  Summary  and  conclusions 

We  have  generalized  previous  analytic  models  for  daylight 
charging  of  a  fast  spinning,  spherical,  dielectric-covered 
spacecraft  to  treat  finite  rotation  rates,  for  the  case  where  the 
sun  direction  lies  in  the  satellite  bellyband  plane.  The  finite 
spin  rate  is  implemented  by  including  an  associated  Legendre 
term  in  the  potential  expansion.  This  term  removes  the  ax¬ 
ial  symmetry  from  the  models  and  allows  differential  charg¬ 
ing  between  back  (shaded)  to  front  (sunlit)  and  also  between 
transverse  surfaces  during  a  rotation.  The  solution  is  devel¬ 
oped  by  introducing  a  rotation  parameter  which  spans  the 
spin  rates  from  the  zero  spin  case  to  the  fast  spin  limit.  Plots 
of  the  solutions  with  varying  spin  show  that  the  rotation  rate 
is  an  important  parameter  in  determining  the  potential  distri¬ 
bution. 

The  models  are  limited  to  daylight  equilibrium  charging 
with  the  sun  direction  in  the  spin  plane.  It  is  assumed  that  the 
ambient  density  is  low  enough  that  Laplace’s  equation  for  po¬ 
tentials  is  valid.  Multipole  terms  higher  than  the  quadrupole 
are  not  considered  in  the  potential  expansion.  Self  consis- 
tant  photoemission  dynamics  is  not  addressed  and  there  is 
a  presumption  of  high-level  charging.  The  satellite  body 
is  assumed  to  be  spherical  and  covered  with  dielectric  film 
of  uniform  resistivity.  The  spacecraft  surfaces  are  not  cou¬ 
pled  via  electrical  connections.  Surface  material  properties 
and  current  collection  algorithms  have  not  been  considered. 
These  many  assumptions  are  sufficient  to  allow  a  solution  for 
barrier  dominated  sunlight  charging  to  be  described  analyti¬ 
cally.  Based  on  the  models,  the  following  (not  all  indepen¬ 
dent)  general  remarks  may  be  made: 

-  The  basic  time  scales  for  charging  are  obtained  by  com¬ 
paring  the  the  dark  side  equilibrium  surface  charging 
time  and  the  satellite  spin  half  period. 

-  “Slow  spin”  occurs  when  the  dark  side  surface  charging 
time  is  short  with  respect  to  the  spin  half  period  and 
“fast  spin”  corresponds  to  a  charging  time  greater  than 
the  half  period. 

-  In  the  potential  expansion,  the  monopole  term  is  always 
present.  The  X-dipole  term  is  largest  at  zero  spin  and 
ramps  down.  The  Y-dipole  term  ramps  up  to  a  maxi¬ 
mum,  then  relaxes  back  down  and  it  occurs  at  all  spin 

www.ann-geophys.net/24/2599/2006/ 


M.  Tautz  and  S.  T.  Lai:  Analytic  models  for  a  spherical  satellite  charging  in  sunlight 


2609 


rates  except  zero  and  infinity.  The  quadrupole  compo-  Appendix  A 

nent  is  zero  at  zero  spin  and  is  the  dominant  term  at  fast 

spin.  Nomenclature 


-  As  the  satellite  spins  up,  the  quadrupole  term  becomes 
stronger  and  brings  in  azimuthal  smoothing  until,  at  fast 
spin,  the  bellyband  potentials  are  all  equal. 


-  The  photoemission  barrier  height  need  only  be  a  few 
volts  to  suppress  photoelectrons,  so  that  for  high-level 
charging  the  solutions  are  close  to  the  threshold  config¬ 
urations.  The  barrier  radius  and  height  tend  to  decrease 
with  increasing  spin. 


-  The  average  surface  potential  is  equal  to  the  monopole 
value.  Changes  to  the  bellyband  potentials  are  linearly 
correlated  with  changes  in  the  potentials  at  the  spin 
poles. 


-  The  differential  voltage  along  the  sun  line  ramps  down 
from  a  maximum  for  a  stationary  satellite  and  goes  to¬ 
wards  zero  when  the  surface  charging  time  approaches 
one  half  the  period  of  a  spinning  spacecraft. 


-  A  differential  voltage  transverse  to  the  sun  line  and  spin 
axis  develops  whenever  there  is  a  departure  from  the 
zero  spin  or  infinite  spin  limit  and  there  is  a  correlated 
lag  angle  between  the  direction  of  the  maximum  sheath 
barrier  and  the  sun  line. 


-  There  is  a  differential  voltage  (azimuth  dependent)  be¬ 
tween  the  spin  poles  and  the  bellyband  surfaces.  In  the 
fast  spin  limit,  it  becomes  azimuthally  symmetric  and 
represents  the  dominant  differential  charging. 


-  At  the  zero  spin  limit  the  monople-dipole  system  is 
recovered  and  at  the  infinite  spin  limit  the  monopole- 
quadrupole  system  occurs.  At  finite  spin  there  is  a 
combination  of  the  monopole,  X  and  Y  dipole,  and 
quadrupole  terms. 


These  are  the  characteristic  features  describing  charging  of 
a  uniform  non-conducting  satellite  spinning  around  the  Z- 
axis  at  a  finite  rate  with  the  sun  direction  along  the  X-axis. 
The  models  describe  schematically  the  shape  of  the  average 
potential  distribution  surrounding  a  charged  spherical  space¬ 
craft  in  a  low  density  space  plasma  and  could  serve  as  a  ref¬ 
erence  point  for  sunlight  charging  of  a  spinning  satellite. 
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relative  coefficients  in  Legendre  ex¬ 
pansion 

model  differential  charging  param¬ 
eters 

transverse  differential  potential 
scaling  function  of  barrier  height 
over  monopole  strength 
radial  distance  from  the  center  of 
the  sphere 

photoemission  barrier  radius  from 
center  of  the  sphere 
photoemission  sheath  width 
polar  angle  function  x=cos & 
function  terms  in  the  an*  coeffi¬ 
cient 

general  coefficients  of  Legendre  ex¬ 
pansion 

spin  rate  dependent  Legendre  coef¬ 
ficients 

sun  aligned  differential  potential 
monopole  potential  strength  param¬ 
eter 

associated  Legendre  polynomial 
dimensionless  spin  rate  parameter 
time  for  equlibrium  surface  charg¬ 
ing 

terms  in  the  coefficients  of  Legen¬ 
dre  expansion 

Cartesian  coordinates  in  the  space 
frame  of  reference 
surface  potential 
photoemission  barrier  potential 
surface  potentials  at  the  belly  band 
(east, west,  back,  and  front  sides) 
surface  potentials  at  the  north  and 
south  spin  poles 

average  potential  at  the  middle 
(bellyband) 

average  spin  pole  potential 
non-monopole  strength  parameter 
polar  angle 
azimuth  angle 
photoemission  sheath  angle 
satellite  spin  period 
potential  limits  for  time  averages 
calculation 
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